%% 脚本说明
%{
- --------------------------------------------------
* 文件名:Contour_main
* 功 能:构造数据，绘制出等高线图
* 变量说明:
-- circle_range 读取文件的名称序列
--
* 注意事项:
可以输入三个变量或者四个变量
需要跳转当次文件的路径才能正常运行，否则请使用绝对路径
- --------------------------------------------------
%}

fclose all; close all;
% var = 'theta';
circle_range = 0:3:12;
% circle_range = [0:3:27,28,29,30];
Input_diretory = "..\Split\Output Results\theta 0-12\09-02-00-17\";
freq_range = [1, 0.397, 0.903]; % 指定频率范围，1表示开，0表示不启用这个功能
YY = circle_range; X = []; Z = []; Y = [];

for eachfile = YY
    myfilename = Input_diretory + num2str(eachfile); %文件名
    tmp = readsweep(myfilename, eachfile);

    if freq_range(1) == 1
        tmp(tmp(:, 1) < freq_range(2), :) = [];
        tmp(tmp(:, 1) > freq_range(3), :) = [];
    end

    X = [X, tmp(:, 1)];
    Z = [Z, tmp(:, 2)];
    Y = [Y, tmp(:, 3)];
end

Paras = {
    'X:', 'Frequency (THz)';
    'Y:', 'Incident Angle \theta (°)';
    'Z:', 'Absorption';
    };
Contour_Plot_YW(Y, X, Z, Paras, [0.0 1]);
